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ABSTRACT 

It is shown that the solution of the quark Dyson-Schwinger equation in 
QCD obtained with a gluon propagator of the form D{q) ~ q'^ /[q'^ + b'^] and 
a quark-gluon vertex that is free of kinematic singularities does not describe 
a confined particle and that there is always a value of b^ = 6^ such that 
chiral symmetry is not dynamically broken for 6^ > b^. 

1. Introduction 

Dynamical Chiral Symmetry Breaking (D^SB) and confinement are two crucial 
features of quantum chromodynamics (QCD). They are responsible for the nature 
of the hadronic spectrum; D^SB ensuring the absence of low mass scalar partners 
of the pion and confinement ensuring the absence of free quarks, for example. A 
natural method for studying both DySB and confinement in QCD is the complex of 
Dyson-Schwinger Equations (DSEs).^ 

One goal of DSE studies is to develop this nonperturbative approach to the point 
where it is as firmly founded as lattice QCD and calculationally competitive. Al- 
though more needs to be done in order to achieve this goal there has been a good 
deal of progress, especially in the study of Abelian gauge theories where direct and 
meaningful comparisons can be made, and agreement obtained, between the results 
of lattice and DSE studi^.i 

Herein a recent studyS of the fermion DSE in which the gluon propagator vanishes 
at = 0, aimed at determining whether such a gluon propagator can support D^SB 
and/or generate a confining quark propagator, is described. There has been renewed 
interest in such a form of the gluon propagator, which was argued in Ref. [^] to be 
associated with the elimination of Gribov copies, because of the recent work of Habel 
et. alB and Zwanziger.El 

In Sec. 2 the DSE for the quark propagator is described in detail. There is also 
a discussion of what is known about the dressed gluon propagator and quark-gluon 
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vertex in QCD. The analysis of D^SB is reported in Sec. 3 and quark confinement is 
discussed in Sec. 4. The results are summarised in Sec. 5. 

2. Dyson-Schwinger Equation for the Fermion Self Energy 

In Minkowski space, with metric g^v = diag(l, —1, —1, —1) and in a general co- 
variant gauge, the inverse of the dressed quark propagator can be written as 

S-\p) =^-m-J:{p) = Z-^{p^)(^^- M{p^)) =A{p^) ^-B{p^), (1) 

with: m the renormalised, explicit chiral symmetry breaking mass (if present); 
the self-energy; M{p'^) = B{jP') / A{p'^) the dynamical quark mass function; and 
Z{p'^) = A'^Ip"^) the momentum-dependent renormalisation of the quark wavefunc- 
tion. The unrenormalised DSE for the inverse propagator is 

S-\p) =^-m^-^-4/| 0^^rS{k)T'^{k,p)D,^{{p-kn (2) 

where D^i,{q'^) is the dressed gluon propagator and is the proper quark-gluon 
vertex, which is illustrated in Fig. |1]. 




Fig. 1. The Dyson-Schwinger equation for the quark self-energy. 
The renormalised, massless (m = 0) DSE can be written as 

J:M = (1 - Zs) ^ + ^Zr\g'j'' 0^,YSn{k)TUk,p)D^AiP ' k)% (3) 

where Zs and Zy are quark-propagator and quark-gluon-vertex renormalisation con- 
stants, respectively, which depend on the renormalisation scale, /x, and ultraviolet 
cutoff, A. Hereafter we suppress the label R. 

The solution of this equation provides information about D^SB . The quark 
condensate, (gg) oc tr[5'(a:; = 0)], is a chiral symmetry order parameter. If there is a 
solution of Eq. (0) with B ^ Q then the quark has generated a mass via interaction 
with its own gluon field and the chiral symmetry is therefore dynamically broken. The 
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solution also provides information about quark confinement, as discussed in Sec. 4. 



2.1a Gluon Propagator 

In a general covariant gauge the dressed gluon propagator, which is diagonal in 
colour space, can be written: 



n(g2 



(4) 



where n(g^) is the gluon vacuum polarisation and ^ is the gauge parameter. In covari- 
ant gauges the longitudinal piece of this propagator is not modified by interactions, 
which follows from the Slavnov- Taylor identities in QCD. 

The Dyson-Schwinger equation for the gluon propagator is given diagrammatically 
in Fig. ^. The symmetrisation factors of 1/2 and 1/6 arise from the usual Feynman 
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Fig. 2. The Dyson-Schwinger equation for the gluon propagator. [Here 
and below the broken line represents the propagator for the ghost field.] 



rules, which also require a negative sign [unshown] to be included for every fermion 
and ghost loop. This equation has been studied extensively.Elo&B There have also been 
attempts to determine the gluon propagator from numerical simulations of lattice- 
QCD.ffl 

The results of the DSE and lattice studies are summarised in Sec. 5.1 of Ref. [|T| 
and are represented in Fig. |^. This figure illustrates that for spacelike-g^ > 1 GeV^ the 
gluon propagator is given by the two-loop, QCD renormalisation group result, with 
the next order correction being <10%. For spacelike-g^ < 1 GeV^, however, the form 
of the propagator is not known. The DSE studies of Refs. suggest a regularised 

infrared singularity, represented by 1 /g^ in the figure. That of Ref. |^ , which differs 
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mainly in that the Ansatz used for the triple-gluon vertex has kinematic singularities, 
suggestsan infrared vanishing form, characterised by / (5^ + 6^), which has also been 
argueduS to be the form necessary to completely eliminate Gribov copies. 

The lattice Landau-gauge studies of Ref. [0 favour the massive vector boson 
form, + m^, which is broadly consistent with the improved studies of Ref. 



On 16^ X 40 and 24^^ x 40 lattices at /? = 6.0 these numerical results allowed a fit of 
the form gVl?" + with h ~ 340 MeV, but a fit using a standard massive particle 
propagator could not be ruled out. On a lattice of dimension 16^ x 24 at /3 = 5.7 
it was found that the gluon propagator was best fit by a standard massive vector 
boson propagator with mass ~ 600 MeV. There is a problem with these studies, 
however, which is indicated by the dashed vertical line at the right of Fig. ^. With 
present technology, the domain of spacelike-g^ < 1 GeV^ is actually inaccessible in 
lattice studies and, since all forms of the propagator are very nearly the same outside 
this domain, it is clear that these results are both qualitatively and quantitatively 
unreliable. 

This brief discussion indicates that at present one can only say that 

2 1 

^ ^ g2 [1 - n(g2)] = TTfe^ 

is not implausible in QCD, at least at small q^. For this reason it is of interest to 
determine whether such a form of D{q^) can lead to D^SB and a confining quark 
propagator, which will provide further insight into the validity of this form of gluon 
propagator. In keeping with the ultraviolet behaviour illustrated in Fig. ^ the results 
obtained with the "ultraviolet-improved" form: 

^D{q^) = a{r-q')^-, (6) 
47r + 0* 

with a(r; q^) = {dir)/ (in [r + q^/AQco]), d = 12/[33 - 2Nf], Nj = 4, where r > 1 
is an infrared regularising parameter, are described herein. 

2.1b Gluon Condensate 

The gluon condensate can be calculated from the nonperturbative part of the 
gluon propagator.0 Using Eq. (0) one obtains 



(a5G'G')^ = -— ln(ln 



A2 



(7) 



For h = 340 MeV, Aqcd = 200MeV and /i = 1 GeV this yields (asGG)^ = 
—0.0047 GeV^, which should be compared with the value inferred from QCD sum 
rules: {agGG) ~ 0.04 GeV^. The "wrong sign" is due to the fact that Eq. (D is 
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weaker than perturbative gluon exchange for all spacelike-g^ and is an harbinger of 
the results to follow. 



2.2 Quark- Gluon Vertex 

The quark-gluon vertex satisfies its own DSE but hitherto there have been no 
studies of this equation in QCD. Making use of the "Abelian approximation",!^ the 
Slavnov- Taylor identity for this vertex reduces to the Ward-Takahashi identity famil- 
iar from QED: 

k,r>^{p,q) = S-\p)~S-\q) , (8) 
k = {p — q), and the quark-gluon vertex can be written in the general formEl 



r,{p, q) = rf (p, g) + ^ r;(p, q) g\p\p- q, q' 



(9) 



1=1 



where the eight tensors, T^, are transverse, {p — q)^T^^{p, Q') = and 

Vf^{p,q) = l:A{v,q)l^. + {p + q)^.[^A{p.q)\b■p + l■q]-^^B{p.q)^m 

with = [F(p2) + F(g2)]/2 andA^fag) = [F{p'^)-F{q^)]/\p^ -q%{oi F = A 

or B. In this approximation QED studieslj can be used to place constraiiits on the 
functions g^. Taking these constraints into account leads to the Ansatz:EJ = 0, 

i 7^ 6, and 



with d{k,p) 



(A;2 + p2 



2d{k,p) 
{e-p'f + 



Aik') - A{p') , (11) 



■ (12) 



This Ansatz, Eq. (D with Eq. (0) and = for z 6, satisfies the Ward- 
Takahashi Identity, is free of kinematic singularities (i.e., has a well defined limit as 
{p ^ q), reduces to the bare vertex in the free field limit in the manner prescribed 
by perturbation theory, transforms correctly under charge conjugation and Lorentz 
transformations and preserves multiplicative renormalisability in the quark DSE. (Of 
these properties the vertex T^'^ satisfies all but the last and, for the most part, the 
calculations described herein were performed with the Ansatz = Fj^*-'.) 

It should be noted that the absence of kinematic singularities is an important 
and physically reasonable constraint. To understand this one can simply consider an 
analogy with PCAC In the case of the axial-current vertex, F^, there is a kinematic 
singularity in the chiral limit, which is identified with the massless pion excitation. A 
kinematic singularity in the quark-gluon vertex would therefore entail the existence 
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of an hitherto unknown, massless excitation in QCD. 

2.3 Quark DSE in Euclidean Space 

In the Abehan approximation one has = 1 = Zr at one-loop in Landau gauge. 
Using this resuh here considerably simplifies Eq. (^), which can be written in Eu- 
clidean space, with metric 5^^ = diag(l, 1, 1, 1) and 7^ hermitian. 



4 , d^k 

l.[p) = 

where 



m = f -0^,i,sik)Uk,p)D,^iip- kY) (13) 



S-^{p) =i-f-p+ S(j9) = i7 ■ pA{p'^) + B{p'^) (14) 

and all the other elements in this equation are taken to be specified by the expressions 
given above evaluated at Euclidean (spacelike) values of their arguments. The im- 
portant, subtle considerations associated with the transformation between Minkowski 
and Euclidean space are discussed in Sec. 2.3 of Ref. [|1|. 

Equation (|l^) is actually a pair of coupled, nonlinear integral equations for A(p^) 
and -B(p^). With = Tj^^ one obtains 

^, 2^ 167r d^k 2n {p-kf 1 ^, 

^(^) = — y (2^"^"' ^^-^^\p- kY + A\k-)k- + B\k-) ^ 

|3i?(fc2)^4(fc2^^A(p2) ^ J^^^2)^^(^2^p2) _ A{e)^B{k\p')\ h{p, A;)| , 

with h{p, k) = 2 [k'^p^ — (A; ■ pY] /{p — k)'^, and a similar but more complicated equa- 
tion for v4(p^), while with = (r|^^ + Tj^ g^), from Eq. ([TID , one has 

B{p^) = RHS of (15) (16) 

levr .A _ {p-kf B{e)^A{ey) {k--p-) _ 

^3 J (27r)4"^ '^^ ^ \p-ky + b^ A^{k^)k^ + B^{k^) 2d{k,p)^ 

with, again, a similar but more complicated equation for v4(p^). 

The results discussed below were obtained from an iterative, numerical solution 
of these equations on a logarithmic grid of x = p^/Aq^p and y = k"^ / K^qj^ points. 
The solutions were independent of the seed-solution and grid choice and also of the 
UV cutoff, which was A^ = 5 x 10* ^\cd- 

3. Dynamical Chiral Symmetry Breaking? 

The gauge-invariant quark condensate is an order parameter for D^SB and is 
obtained from the trace of the quark propagator: 
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where is the renormahsation point for the condensate, which is usually fixed at 
1 GeV. This is the parameter that is used to study DxSB in lattice QCD. A nonzero 
value signals D^SB. 

In Fig. ^ the condensate obtained from the numerical solution of Eq. (|15D and 
the associated equation for A(p^) for values of Inr in the domain [0.0, 0.7] and 6^ in 
[0.1,1.0] is plotted. This figure shows regions of unbroken and dynamically broken 
chiral symmetry. Notably there is no D^SB for the value of 6^ ~ 3 Aq(jj3 inferred 
from lattice simulations. 

The numerical results suggest that {qq)ij. rises continuously from the transition 
boundary and hence that the transition is second order. Assuming therefore that, for 
a given value of Inr, the order parameter behaves as ~ C (1 — b'^/bl)^ one 

obtains (3 = 0.572 with = 0.020. Including the term in the vertex only leads 
to a small quantitative change in the results. For example, with Inr = 0.6 one finds 
(3x6 = 0.579, cT/3^g = 0.015. This suggests that for any vertex that is free of kinematic 
singularities there is a critical value of 6^ = bl such that there is no D^SB for b'^ > b^. 



4. Quark Confinement? 

In order to determine whether the quark propagator obtained as a solution to 
Eq. ([T3|) with an infrared- vanishing gluon propagator can represent a confined particle 
we follow Ref. |l^ and adapt a method commonly used in lattice QCD to estimate 



bound state masses. Writing <Js{s) = B{s)/[sA{sY + B{sY], s = p^, defining 



and, for notational convenience, E{T) = — ln[As(T)], it follows from the axioms of 
field theory that if there is a stable asymptotic state with the quantum number of the 
quark then 

dE(T) , , 

lim ,^ ' = m ; 19 

where m > is the mass of this excitation; i.e, this limit yields the dynamically 
generated quark mass. 

As a simple example one can consider the Nambu-Jona-Lasinio model0 in which 
the dressed quark propagator is S{p) = ■ p + M] and hence (Js{s) = M/[s + M^]. 
In this case As(T) = exp(— MT) /2, which, from Eq. (|I9|) , yields m = M, as one 
would expect. 

This confinement test was applied to the numerical solutions in the following 
cases: 1) The propagators obtained with Inr = 0.1 and 6^ in the range [0.1, 1.0]; 2) 
The propagator obtained with Inr = and 6^ = 0.35, which yields the largest value 
of {—{qq)^) on the (6^, Inr) domain considered; 3) Two propagators obtained with 
(62, Inr) = (0.1,0.6) - one using = Fj^^ and the other using F^ = (Fj^^ + /). 
Plots of E'{T) for the family of propagators obtained with Inr = 0.1 are presented in 
Fig. ^ (It should be noted that a quark propagator with complex conjugate polesi 
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leads to a form of As(T) with zeros and hence to E'{T) with zeros and poles; a 
strong signal of which there is no sign in Fig. ^) Since the behaviour of all the other 
solutions was qualitatively the same as that described by the results presented in 
Fig. H one concludes that, independent of b, D{q) ~ q'^/lq'^ + b'^] does not yield a 
confining quark propagator. 

5. Summary 

A study of the quark Dyson-Schwinger equation in QCD using a gluon propagator 
that vanishes as — > 0, Eq. (|), and a dressed quark-gluon vertex, T^{p, q), that has 
no kinematic singularities; i.e., has a well defined limit as p —>■ q, is described. The 
results indicate that such a gluon propagator can only support D^SB for values of 
Inr and 6^ less than certain critical values and does not confine quarks. The results 
are quahtatively independent of the model dependent elements in this study. One is 
therefore lead to conclude that the dressed gluon propagator in QCD does not vanish 
in the infrared. 
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Fig. 3. Results of studies of the gluon propagator in QCD. Typical 
values of nig — 0.5 GeV and b — 0.7 GeV have been used. 
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Criticality plot, -<qq>^^ vs. b^, ln(T) 




0.3 
b2 



Fig. 4. Criticality plot for (—{qq)^)^^ as a function of Inr and P. The 
condensate, {—{qq)^)^, is in units of MeV, scaled to fi^ — IGeV, and 
is in units A.qq^ [ — 0.49 =^ 6 ~ 140 MeV]; the gluon regulator r is 
dimensionless. 
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Mass curves, ln('i:)=0.1, Minimal Ball-Chiu Vertex 




Euclidean time T 

Fig. 5. Dressed-quark-mass curves for the family of propagators with 
the minimal Ball-Chiu vertex and Inr = 0.1. The masses are in units of 
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